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CONFIDENTIAL BULLETIN 


INTERNAL AND EXTERNAL AERODYNAMICS OF DUCTED 
BODIES AT SUPERSONIC SPEEDS 


By Clinton E. Brown 


SUMMARY 


A method for calculating the lift, moment, and 
pressure drag of slender open-nose bodies of revolution 
at supersonic speeds is described. An application of the 
method to a typicel ram- jet fuselage is shown to give 
excellent agreement with available experimental data. 

A drag comparison was omitted because of the presence 
of skin-friction drag In the experimental drag data. 

The problem of obtaining high tot al-pre3sure reoovery at 
supersonic speeds Is discussed and some experimental 
data obtained at the Langley Memorial Aeronautical 
Laboratory of the NACA on circular diffusers Is pre- 
sented. It is oolnted out that variable-geometry dif- 
fusers might be designed to give high diffuser efficien- 
cies over a wide range of Mach number. 


INTRODUCTION 


Tflfork on ram- jet housing bodies was started at the 
NACA In the summer of 19^3* At that time, no experi- 
mental data of any kind were available in this country 
on ducted bodies at supersonic speeds. As work pro- 
gressed, an open-nose body of revolution was selected 
as a fuselage for a possible supersonic airplane. A 
test f.rogram was started in the Langley model super- 
sonic tunnel and a 5~inch model of the fuselage selected 
was tested at several Mach numbers. Curves from these 
data were obtained of lift, moment, and drag coefficients 
In oonjunotior with the test program, an investigation 
was begun to find means for calculating the aerodynamic 
characteristics of these bodies. As a result, a method 
was found (reference 1) which, for the body tested, 
gave excellent agreement with the experimental results. 
Previous to this work, a series of tests were conducted 
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in a 1-Inoh supersonic jet at the Langley Memorial 
Aeronautical Laboratory •on-oiroular diffusers having a 
small contraction Just behind the duct entrance (refer- 
ence 2). Some of - the results obtained are presented 
and discussed. 

The discussion given in the present paper was 
originally presented at a Symposium on Supersonic Plow 
held at the Johns Hopkins University, ADplied Physics 
Laboratory, Silver Spring, Md., on Dec. S and 7» 19^5* 


SYMBOLS 


x, r, 0 cylindrical coordinates 

X distance along X-axis measured from nose of 

body 

l length of body 

R radius of body 


P 


Mach angle 



B = v^ 2 “ 1 
0(x, r, 0) 
#l(x, r) 
0 2 (x, r, 0 ) 


perturbation potential 
perturbation potential for axial flow 
perturbation potential for cross flow 


axial velocity increment 
radial velocity increment 



V 


velocity In undisturbed stream 
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a 

M 

P 

A? 

P 

a 

0 



C 


D 


C 


m 


velocity of sound in undisturbed stream 

Mach number -in undisturbed stream (V/a) - - 

density in undisturbed stream 

incremental surface pressure due to angle of 
attack 

local pressure 

pressure in undisturbed stream 

ratio of soecific beats of air (l.U) 

angle of attack, radians (except where otherwise 
noted) 

angle between surface of body and X-axis 


lift coefficient / ■ 
drag coefficient f 


1) 


moment coefficient 


( Drag \ 

tw) 

/ foment \ 


u variable of integration 
<| = x - Br cosh u 

Jl = Xf - B p i 



*n ~ Jj 
BR n 


Ai = f 1 (S)i 


k 
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Subscripts and superscripts:' 

< . 

N refers to nose 

n refers to n tla integration station, summation 

variable 

i refers to i*'* 1 integration station, summation’ 

variable 

max maximum 

deg in degrees 

MATHEMATICAL* ANALYSTS 
Pointed 3odies of Revolution 


The analysis that is presented herein is essentially 
that found in reference 1. The mathematical methods used 
were first developed by von Harman arid Moore (refer- 
ence aoplied to the cone problem and to sharp* 

nose projectiles. If the problem is restricted to 
slender bodies, the differential equation of motion for 
a compressible fluid can be linearized to give. In 
cylindrical coordinates, the simple form 


dr r dr r 2 £>9 . d* 


where /f is the potential function assumed to represent 
the effect of a small disturbance set up by the slender 
bodies being considered. The problem is, therefore, to 
find a solution of equation (1) that will satisfy the 
known boundary conditions at the surface of the body. 

A general solution of the differential equation (1) 
when M > 1 for diverging waves has been found by Lamb 
(reference *;.) to be, with a slight change in notation: 


s 


ftgr s cos s8 + P s r s sin 89 



* 
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whore 



p s = (db/ p o 

> 


and 





09 

f(x - Br cosh u) du 


P 


0 



Br cosh u) du 


where 


B = 



The part of Lamb's general solution corresoondl ng to 
converging waves does not anoly to tho presciio oroblem 
since all disturbances originate on the body av ’ diverge 
Into the flow field Investigated. Von Farman and Koore 
have investigated the problem cf the resistance of pro- 
jectiles and cones (reference 3) and have found a 
solution for the case of axial symmetry 

$ - - \ f ( x - Br cosh u) du (3) 

Jo 

which oan be seen to be a special case of the general 
solution with s = 0. In their analysis It was found that 
the body, in this case a sharp-nose projectile, could be 
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represented by q distribution of sources along the s ' 

stax}£|gg‘a£ the nose pf the body. By a numerical 

method pf £fitegF«t£iqn £he equations for the velocity 

increments y ‘ and v could then be written 

j x 


n = - B ^ A 1 y(*l °) Z - 1 - /( T l-I n ) 2 - X j 


( 4 ) 


*n 


= 51 A i ("oosh- 1 ^ 11 ) - cosh.- 1 (t 1 . 1 “) 
1=1 L 


( 5 ) 


where 

T n _ *n ~ 

1 

and 

ii = *i - BR i 


A i = r ’ <£>* 


with the boundary conditions 


V + v 


dr 

dx 


( 6 ) 


These three 


equations in three unknowns 



at each station on the body for v„ and v„. 

" r x 


were solved 
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The pressures were then found from the Bernoulli equation 
in the form: 


_ 1 _ 



Ferrari (reference 5 ) and Tsien (reference £) have 
Independently found solutions for the case of nointed 
bodies of revolution at small angles of attack. It was 
shown that the potential could be expressed in two terms: 
the first, from equation ( 3 )* 


r 

0i ~ ~ I r i- x - Br cosh u) du 

Jq 

is the solution for the cure axial flow already described, 
and the second 



-3 cos 



Br cosh u) cosh u du 


( 8 ) 


represents the cross-flow potential of an arbitrary 
distribution of doublets along the axis of the body 
starting at the nose of the oone or projectile. The 
form of equation ( 6 ) is such that the cross flow is 
from the direction 0 = 0 . 

By neglecting the small effect of the axial flow 
on the lifting pressures, Tsien obtained for the pointed 
projectile of arbitrary shape the equations: 
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The values in these equations are assumed to he 

constants for each interval of the step-by-step process. 
The moment coefficient of equation (10) is assumed posi- 
tive for nosing-up moments, the moments being taken 
about the nose. 
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Open-Hose Bodies 

The flow -conditions over an open-nose body differ 
from those of pointed bodies In that, for finite angles 
of the. nose lip, the flow Is two-dimensional at the lip* 
Oils problem was not considered In references 3, 5, 
and o, and the general solution should therefore be 
examined, to determine Its applicability to this speolal 
case. It has been shown by Lamb (reference Ij.) that a 
sufficient requirement for the existence of the general 
solution to the differential equation of motion is that 
f(x - Br cosh u) be zero for all values of the argument 
less than some arbitrary limiting value* Bio determina- 
tion of f (x - Br cosh u) such that the boundary condi- 
tions at the open-nose body are satisfied assures the 
fulfillment of this general requirement. For the usual 
case of supersonic flow into the nose, the boundary 
condition requires the surfaoe of the body to be a 
continuation of a cylindrical stream surfaoe of radius Rjf 
in the undisturbed flow ahead of the body, as shown in 
figure 1. The perturbation potentials (equations ( 3 ) 
and ( 8 )) must therefore be zero at the cylindrical stream 
surfaoe ahead of the body. Substituting § * x - Br cosh u 
in equations ( 3 ) and ( 8 ) gives 




y/( x — g ) 2 - E^r 2 


( 12 ) 


and 



00 s 3 


px-Hr 


f 2 (£)(x - £)d£ 

>/ (x - - B^r^ 


(13) 


The boundary conditions are obviously satisfied when 
f^(£) = fg(£) = 0 for all values of £ < Xq - BRg, 

where the point (xq, Rq) Is at the lip of the open-nose 

body. Values of f-j^S) and f 2 (£) fox* £ > Xq - BR^ 
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then remain to be determined so that the bod; surface 
Is a continuation of this stream surface. From physioal 
considerations, 1 f^CS) hnd f 2 (S) may be regarded as 

an social distribution of sources and doublets, respec- 
tively, where £ is measured along the X-axis. As the 
effect of a source or doublet can only be felt along or 
behind its Jfach cone, the source distribution must begin 
a distance BB Q ahead of the, nose- This point is 

chosen for the origin of the coordinate" system, it must 
be emphasized that the source and doublet' distribution 
determined by satisfying the boundary conditions at the 
stream and body surfaces shown in" figure 1 does not 
represent oorrectly the flow infside that stream surface. 
This result corresponds to the physical fact that the . 
actual supersonic flow into the nose does 'not affect 
the flow external to the body.. The basic assumptions 
of potential flow and small disturbances are valid 
provided the slope of the body surface Is small. Actually, 
for finite angles of the nose lip, a nonconioal shock 
wave is formed which causes a loss In total head and 
produoes rotation in the field. . 

Numerical integration of equation. ( 12) , with con-r 
stant values of ^'(S) over the integration Intervals 

assumed, results In the same expressions for v and v x 

r n n 

as those obtained by von Kdrmdn r and Jioore ( equations ( 4 ) 
and ( 5 ))» These-censtant values of f^!(§) are deter- 
mined by satisfying the boundary condition: 


— — = tan 6 (111.) 

V + n 


Triers tan 6„ is the slooe of surface of the body at 

XiV II 

the n txl interval of integration. 

By following the method of reference 6, the lift 
and moment coefficients for small angles' of attack 
based on the area of the nose may be written 
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^ nl+BRjj; 


C L = 


2' ■ 


boa 0 <30 R dx 


'0 U BR, 


. I+SRH ( . ^ 
» 2 » I L v s n;6x 


cos 0 d0 R dx ( l6) 


'0 «BR 


where Z is the length of the body, R^ is the nose 

radius, and the moments are taken about the midpoint of 
the body. Substituting the expression for d02/& x - , - n 

equations ( 15 ) and (l6) gives, for C L and G^, 


_ 1>.b 


tg_ f 2 n „„ r l+E \ „ 

*-?- I oob^B d 0 / ■ R dx 
F V 4 ^BR._ 


k x-BR 


f 2 '(Z)(x ~ g)dg 

3E \f( x - l) 2 - E^P 2 


(17) 


C m = 




coa^Q d0 


pl+BR N 

I («-*-») 

v no 


R dx 


,x-BR 


- g)dg 

BR \/(* - g) 2 - B 2 R 2 


(18) 


The distribution function fp(£) must be determined by 
the boundary condition ^ 
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(19) 


which assumes that the radial velocity is normal to the 
surface. A more rigorous boundary condition taking 
into account the slope of the body was given by Ferrari 
(reference 5)* For small surface angles, however, equa- 
tion (19) is within the accuracy of- the small-nerturbation 
assumptions. The expression 


/ dfffgN _ eos 9 P 

vw r=R - r2 / 


cos 0 / ,X " ER f£)f* - £) 2 d£ 


(20) 


\J( x - &) 2 - B 2 !? 2 


is integrated numerically for constant values of 
fg'CS) = % over the 1™ interval of integratio: 


on to 


obtain the sum 


(' 


6 02 ^ _ ^cos 9 

“2 

i=l 


br ) 

' r=H 


y~ N 1 cosh“ 1 (T 1-1 n ) - cosh' 1 ^ 11 ) 


+ ( T i-i n ) \/( T i-i n ) 2 - 1 - T i n W^f- 

Substituting this equation in equation (19) gives 

1 = t. [»»-*“ ■ 1 ( T 1 -i n ) - oosh ' 1 (k n ) 


(21) 


+ (Ti-i n )\/(Ti-i n ) - 1 “ T i n \/( T 3 n ) - 1 


( 22 ) 
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In equations ( 23 ) and ( 2 k) the pressure used for a 
given integration interval is the average of the pres- 
sures at the beginning and at the end of the interval. 
This scheme' of using average lifting pressure is particu- 
larly necessary in regions where the pressure is rapidly 
changing. The method does not give the pressure at the 
beginning of the first integration interval, that is, 
at the point n = 0. It can be shown that, as the first 
interval approaches zero, the pressure at ths .lip (n = 0 ) 
is obtained by letting the expression in equations ( 23 ) 
and ( 2 k) 



have the value 0.5 when n = 1 . 


COMPARISON OP THEORY AND EXPERIMENT 


Calculations were made to obtain the pressure drag, 
lift, and moment of the open-nose body on which tests 
had been made. Figure 1 shows a diagram of the body 
for which the calculations were made, with the integra- 
tion stations at which the boundary conditions were 
applied. The computations were made for the Mach 
numbers l.k 5 » 1 » 66 , 1 . 75 * and 3 » 00 . 

The experimental drag results had a considerable 
amount of skin drag and a drag comparison was not there- 
fore conclusive. Estimates of the skin-friction coeffi- 
cient showed, however, that the calculated drag was of 
the right order. Eie agreement of the experimental lift 
and moment coefficients with the calculated restilts was 
quite good. Figure 2 gives a comparison of the lift- 
and moment-ooeff icient curves with the experimental 
data. The effect of the Internal flow on the lift can 
he seen from momentum considerations to be AC^ = 2 a, 

whereas the moment which is taken about the midpoint of 
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the fuselage Is AC m = a. These increments have been 
added to the calculated results for the comparison. In 
the tests, a tail-surface supporting ferrule was placed 
over the' tail section - ;' and" no experimental differences 
in lift and moment could be detected. It was concluded, 
therefore, that the thickening of the boundary layer 
ahead of the sho-ck at the trailing edge was (in effect) 
forming a ferrule tail when the ferrule was not actually 
In place. Additional calculations were then made for 
the shape tfith the ferrule tail and, as can be seen In 
figure 2, the agreement with experiment was Improved. 

The general results of the computations are inter- 
esting as they show the effect of Mach number on the 
body characteri stlos. In figure 3 the lift-curve slope 

Gl / a deg a ^ ows an Increase with increasing tf.ach number 

as Is the case with projectiles, tests of which were 
made in Germany and Italy and reported in reference 7* 
The pressure-drag coefficient C D drops off, as would 

be expected. The effect of -Jach number on moment and 
lift can be better illustrated by plotting the effective 
lifting pressures for different Mach numbers. In fig- 
ure ij. the distribution of Incremental surface pressure 
coefficients 008 0 I s plotted for M = 1.6 

and M = 3*°« As tne Mach number Increases, a greater 
proportion of the lift is shown to be carried over the 
center section of the body and, as a result, a rearward 
shift occurs In the center of pressure. Figure 5 shows 
the distribution of drag pressures over the body. It 
can be seen that over the conical nose section the pres 
sures fall from the lip-wedge pressures (pressures on a 
two-dimensional wedge of the same angle, reference C) to 
approximately the pressure expected on a cone of the 
same apex angle (reference ~j J. At the corners the pres- 
sures fall approximately In accordance with the Prandtl- 
Meyer relation for flow around a corner (reference J) . 
The method of calculation oresented herein is given In 
more detail in reference 1. 

Ferrari (reference 10) has developed the method of 
characteristics for determining the fields of flow about 
arbitrary bodies of revolution. Bis method consists of 
a step-by-step determination of .the pressures and ■veloci 
ties at all points In the field starting from some known 
boundary. The work Involved In obtaining lift and 
moments is, however, especially tedious as the step-by- 
step process must be carried out along and around the 
body. 
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DIFFUSERS 


The problems of obtaining high ram recovery in 
supersonic flows is of great importance. At high Mach 
numbers the pressure recovered by a normal shock wave 
becomes a small fraction of the total pressure in the 
stream; the total pressure losses through shock waves 
can be reduced only by a reduction in the Mach number 
normal to the wave. Interest is therefore centered in 
the design of diffusers in which the intensity of the 
■shock waves is kept to a minimum. Theoretically a 
diffuser could be so shaoed as to allow smooth isen- 
tropic compression through the speed of sound; however, 
F.antrowitz and Donaldson (reference 2) have shown that 
such a flow is unstable and unattainable in practice. 

If, however, a reversed Laval nozzle in which such a 
flow exists were considered, the Mach number at the 
minimum-area soction would then be unity, corresponding 
to a maximum mass flow per unit area. Any disturbance 
in the stream causing a gain in entropy (loss in total 
head) ahead of the minimum will cause a decrease in the 
quantity pV rfax . This decrease in mass flow at the 

minimum leads to an accumulation of fluid which will 
Immediately cause a normal shock to progress upstream 
and the mass flow Into the diffuser will .thus be reduced. 
If the area of the minimum section Is Increased, the 
mass flow can be increased and the shock wave out in 
front will move down to a position In the diverging 
passage. 

Circular diffuser s.- The work of Kantrowitz and 
Donaldson (reference &T was done on circular diffusers 
for the purpose of obtaining designs which, for a 
specified Mach number, would allow the shock to occur 
In a region of reduced Mach number, thereby reducing 
the total-pressure losses through the shock and giving 
more efficient diffusion. The results of the tests 
showed that by a suitable choice of contraction ratio 
(ratio of entrance area to minimum area) the efficiency, 
even wi th a long subsonic diffuser, could be made to 
equal and to exceed, at high Mach numbers, the efficiency 
obtained by a normal shock alone. Figure 6 shows the 
test results on three diffusers of different contraction 
ratios. The efficiency shown m this figure is defined 
as the ratio of the kinetic energy that the diffused air 
would have after Isentropic expansion to free-stream 
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pressure to the kinetic energy of the free stream before 
entering the diffuser. The curves for the diffusers 
that are designed for the lower' 14aoh numbers fall below 
the normal shock lirle. This result, however, is prabably 
the effect of the larger ratio of friction loss to shock 
loss in the lower supersonic Mach number range. The 
curve labeled "Grocco’s r^sUm^" is one given in refer- 
ence 11 as the probable efficiency of diffusers for 
supersonic wind tunnels that have large, well-established 
boundary layers ahead of the shock. 

Annular Inlets .- Several proposals have been made 
for missiles with annular inlets located back on a 
pro Jectile-shaped nose. Inlets of this type seem to have 
several outstanding disadvantages. The most serious is 
the large boundary layer which, in passing through the 
shock wave in the diffuser, will most certainly cause 
serious separation. Another problem is that the inlet 
may be located in a high local Mach number region in 
which the diffuser efficiency would be again impaired. 

A preliminary test of such a diffuser has been made in 
the Langley model suoersonic tunnel. As the tests were 
cut short, only a few points were obtained; however, the 
best efficiencies found were consideraoly below those 
obtained by Kantrowltz and Donaldson. Total- he ad 
recovery of the model was 'JO oercent at a ^ach number 
of 1.55* In the work of Taylor and fc.accoll (refer- 
ence 12) on flow over cones, it was found that for large 
cone angles subsonic flow was produced on the surface of 
the cone. This result immediately suggests the possibility 
of an annular diffuser built around a cone. At present, 
further work is in progress at the Langley Laboratory of 
the NACA on diffusers of this type. Recent reports of 
work done In Germany on similar diffusers have indicated 
that substantial gains in efficiency have been oDtained, 
especially at high Mach numbers. 

Variable-area diffusers .- The use of constant- 
geometry diffusers operating at other than design con- 
ditions will probably always result in reduced diffuser 
efficiencies. If a variable minimum is provided in the 
diffuser, the local Mach number at which the shock wave 
occurs can be controlled and kept as close to 1 as 
stability conditions will permit. It appears that the 
pressure fluctuations In the combustion chamber of Jet- 
propelled supersonic, aircraft will probably limit the 
shock Mach number to some value slightly above 1 where 
the shock losses are still insignificant. The exact 
value of this local Mach number must be determined by 
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tests of the burners and air components In the ducts 
behind the shock wave. A movable cone In, a fixed 
cylinder similar to the area-changing device used on 
several jet-motor exits appears to offer possibilities 
for use as a variable-geometry diffuser; it is also 
possible that a rubber area-changing device similar in 
operation to wing de-icer boots would make a simple 
easily controlled variable diffuser. 


CONCLUDING REMARKS 


Flight in the supersonic speed range will undoubtedly 
require airplanes of high aerodynamic refinement. At 
present, a large and increasing amount of theoretical 
work is available on this subject. The lift and moment 
characteristics of ducted bodies rt supersonic speeds, 
obtained by the method presented, were shown to give 
excellent agreement with available experimental data. 

A great need is evident, however, for more test data at 
high Mach numbers, since theory in this field Is well in 
advance of experiment . 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley field, Va. 
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Figure 1.- Location of integration stations and intervals on ram-jet body. 
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